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Abstract 

Context. The carbon-rich asymptotic giant branch star IRC +10216 undergoes strong mass loss, and quasi-periodic enhancements 
of the density of the circumstellar matter have previously been reported. The star's circumstellar environment is a well-studied, and 
complex astrochemical laboratory, with many molecular species proved to be present. CO is ubiquitous in the circumstellar envelope, 
while emission from the ethynyl (C2H) radical is detected in a spatially confined shell around IRC +10216. As reported in this article, 
we recently detected unexpectedly strong emission from the N = 4 - 3, 6 - 5, 7 - 6, 8 - 7, and 9-8 transitions of C 2 H with the 
IRAM 30 m telescope and with HerschelfHIFl, challenging the available chemical and physical models. 

Aims. We aim to constrain the physical properties of the circumstellar envelope of IRC +10216, including the effect of episodic mass 
loss on the observed emission lines. In particular, we aim to determine the excitation region and conditions of C2H, in order to explain 
the recent detections, and to reconcile these with interferometric maps of the N = 1 - transition of C2H. 

Methods. Using radiative-transfer modelling, we provide a physical description of the circumstellar envelope of IRC +10216, con- 
strained by the spectral-energy distribution and a sample of 20 high-resolution and 29 low-resolution CO lines — to date, the largest 
modelled range of CO lines towards an evolved star. We further present the most detailed radiative-transfer analysis of C2H that has 
been done so far. 

Results. Assuming a distance of 150pc to IRC +10216, the spectral-energy distribution is modelled with a stellar luminosity of 
1 1300 Lq and a dust-mass-loss rate of 4.0 x 10~ 8 Mq yr~' . Based on the analysis of the 20 high-frequency-resolution CO observa- 
tions, an average gas-mass-loss rate for the last 1000 years of 1.5 x 10~ 5 Mq yr~' is derived. This results in a gas-to-dust-mass ratio of 
375, typical for this type of star. The kinetic temperature throughout the circumstellar envelope is characterised by three powerlaws: 
Tten(r) oc r~ - 58 for radii r < 9 stellar radii, 7\j n (r) oc r -° A0 for radii 9 < r < 65 stellar radii, and 7\j„(r) oc r ~ L20 for radii r > 65 stellar 
radii. This model successfully describes all 49 observed CO lines. We also show the effect of density enhancements in the wind of 
IRC + 10216 on the C 2 H-abundance profile, and the close agreement we find of the model predictions with interferometric maps of 
the C2H N = 1 - transition and with the rotational lines observed with the IRAM 30 m telescope and HerschelfHWl. We report on 
the importance of radiative pumping to the vibrationally excited levels of C2H and the significant effect this pumping mechanism has 
on the excitation of all levels of the C 2 H-molecule. 

Key words. Stars: individual: IRC + 10216 — stars: mass loss — stars: carbon — astrochemistry — stars: AGB and post-AGB — 
radiative transfer 



1. Introduction 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



The carbon-rich Mira-type star IRC +10216 (CW Leo) is located 
at the tip of the asymptotic giant branch (A GB), where it loses 



mass at a high rate (~1 — 4 x 10 3 yr 
1 19971: iGroenewegen et all 119981: ICernicharo et al 



Crosas & Menten, 



2000; 



1 
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iDe Beck et all 201Q[). Located at a distance of 120 - 250pc 



(Loud 

GUI" 



(lLoupet al.Lll993blC rosas & MentenL[T997HGroenewegen et all 



Cernicharo et all 120001) . it is the most nearby C-type 
AGB star. Additionaly, since its very dense circumstellar 
envelope (CSE) harbours a rich molecular chemistry, it has 
been deemed the prime carbon-rich AGB astrochemical 
laboratory. More th an 70 molecular speci es have already 



been detected (e . g. ICernicharo etail l2000t iHe et all 120081 
iTenenbaum et all 1201 oh . of which many are carbon chains, 
e.g. cyanopolyynes HC„N (n=l, 3, 5, 7, 9, 11) and C„N (n=l, 
3, 5). Furthermo re, several anions have been identified, e.g . 



Cernicharo eTalT. 120071: iRemii an et all 120071: 
20071). and C yN (iThaddeus et all 120081). 

dAgundez et all 



and C ^N 
dCernicharo et all 120081) . and CN 



C„H (n = 4,6,8 ; 
iKawaguchi et all 

C 5 N- 

l2oToh . Detections towards IRC +10216 of acetylenic chain 
radicals (C„H), for n=2 up to »=8, have been report ed by e.g. 
iGuelin et all Jl978l CdH). ICernicharoetal] d!986alfl Il987bl 
CsH). ICernicharo et all (fl987al) a nd IGuelin et all dl987i C ft H) 
IGuelin et al] (I1997L CtH). and ICernicharo & Guelinl (1 19961 
C 8 H). 

The smallest C„H radical , ethyn yl (CCH, or C2H), was 
first detected by iTucker et all d 19741) in its N = 1-0 tran- 
sition in the interstellar medium (ISM) and in the envelope 
around IRC +10216. It was shown to be one of the most abun- 
dant ISM molecules. The formation of C2H in the envelope of 
IRC +10216 is attributed mainly to photodissociation of C2H2 
(acetylene) , one of the most pro minent molecules in carbon-rich 
AGB stars. iFonfria et al.l d2008l) modelled C2H2 emission in the 
mid-infrared, which samples the dust-formation region in the in- 
ner CSE. The lack of a permanent dipole moment in the lin- 
early symmetric C2H2-molecule implies the absence of pure ro- 
tational transitions which typically trace the outer parts of the 
envelope. C2H, on the other hand, has prominent rotational lines 
that probe the chemical and physical conditions linked to C2H2 
in these cold outer layers of the CSE. It has been established that 
C2H emission arises from a shell of ra dicals situated at ~15" 
from the central star (IGuelin et al.Lll993h . Observations with the 
IRAM 30 m telescope and Herschel/HIFI show strong emission 
in several high-Af rotational transitions of C2H, something which 
is unexpected and challenges our understanding of this molecule 
in IRC +10216. 

We present and discuss the high-sensitivity, high-resolution 
data obtained with the instruments on board Herschel 
dPilbratt et all 1201 Ol): HIFI (Heterod yne Instrument for the 



Far Infrared; Ide Graauw eTatl l2010l), SPIRE (S pectral and 
Photometric Imaging Receiver; Griffi n et all 1201 Ol). and PACS 
(Phot odetector Array Camera & SpectrometenTPoglitsc h et all 
120101) and with the IRAM 30 m TelescopdUin Sect. [2] The phys- 
ical model for IRC +10216's circumstellar envelope is presented 
and discussed in Section [3] The C2H molecule, and our treat- 
ment of it is described in Section [4] A summary of our findings 
is provided in Secfion[5] 

2. Observations 

2.1. Herschel/HIFI: Spectral Scan 

The HIFI data of IRC +102 16 presented in this paper are part 
of a spectral line survey carried out with HIFI's wideband spec- 
trometer (WBS; Ide Graauw et alll2010l) in May 2010, on 3 con- 
secutive operational days (ODs) of the Herschel mission. The 



scan covers the frequency ranges 480 - 1250 GHz and 1410 - 
1910 GHz, with a spectral resolution of 1 . 1 MHz. All spectra 
were measured in dual beam switch (DBS; Ide Graauw et all 
1201 01) mode with a 3' chop throw. This technique allows one to 
correct for any off-source signals in the spectrum, and to obtain 
a stable baseline. 

Since HIFI is a double sideband (DSB) heterodyne instru- 
ment, the measured spectra contain lines pertaining to both 
the upper and the lower sideband. The observation and data- 
reduction strategies disentangle these sidebands with very high 
accuracy, producing a final single-sideband (SSB) spectrum 
without ripples, ghost features, and any other instrumental ef- 
fects, covering the spectral ranges mentioned above. 

The two orthogonal receivers of HIFI (horizontal H, and ver- 
tical V) were used simultaneously to acquire data for the whole 
spectral scan, except for band 2eq Since we do not aim to study 
polarisation of the emission, we averaged the spectra from both 
polarisations, reducing the noise in the final product. This ap- 
proach is justified since no significant differences between the H 
and V spectra are seen for the lines under study. 

A detailed description of the obser vations, and of the data re- 
duct ion of this large survey is given bv ICernicharo et all d2010bl) 
and ICernicharo et al.l (1201 lbl in prep.). All C2H data in this 
paper are presented in the antenna temperature (7^) scal^B 
iRoelfsema et al.l(l2012l submitted) describe the calibration of the 
instrument and mention uncertainties in the intensity of the order 
of 10%. 

For the current study we concentrate on observations of CO 
and C2H. The ten CO transitions covered in the survey range 
from J = 5 - 4 up to / = 1 1 - 10 and from J = 14 - 13 up to 
7=16-15. The C2H rotational transitions N = 6 - 5, 7 - 6, and 
8 - 7 are covered by the surveys in bands la (480 - 560 GHz), 
lb (560 - 640 GHz), and 2a (640 - 720 GHz), respectively. The 
N = 9 - 8 transition is detected in band 2b (720 - 800 GHz) of 
the survey, but with very low S/N; higher-A^ transitions were not 
detected. A summary of the presented HIFI data of C2H is given 
in TableQ] the spectra are shown in Fig.Q] 

2.2. Herschel/SPIRE 

In the framework of the Herschel guaranteed time key 
programme "Mas s loss of Evolved StarS" (MESS; 
iGroenewegen et ail 1201 ll) the S PIRE Fourier Transform 
Spectrometer (FTS; iGriffin et all 1201 Ol) was used to obtain 
IRC+10216's spectrum on 19 November 2009 (OD 189). 
The SPIRE FTS measures the Fourier transform of the source 
spectrum across two wavelength bands simultaneously. The 
short wavelength band (SSW) covers the range 194 - 313/iin, 
while the long wavelength band (SLW) covers the range 
303 - 671 pm. The total spectrum covers the 446 - 1575 GHz 
frequency range, with a final spectral resolution of 2.1GHz. 
The quality of the acquired data permits the detection of lines as 
weak as 1 - 2 Jy. For the technical background and a description 
of the data-redu ct ion pro cess w e refer to th e discus sions by 
ICernicharo et al.l d2010al) and iDecin et al ] d2010bl) . These 
authors report uncertainties on the SPIRE FTS absolute fluxes 
of the order of 15 - 20% for SSW data, 20 - 30% for SLW data 
below 500 pm, and up to 50% for SLW data beyond 500 pm. 



1 Based on observations carried out with the IRAM 30 m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



2 Because of an incompleteness, only the horizontal receiver was 
used during observations in band 2a. Supplementary observations will 
be executed later in the mission. 

3 The intensity in main-beam temperature T MB is obtained via T MB = 
T^/tjmb, where 77MB is the main-beam efficiency, listed in TablefTJ 
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Figure 1. C 2 H N = 1 - 0, 2 - 1, 3 - 2, 4 - 3, 6 - 5, 7 - 6, 8 - 7, and 9 - 8 rotational transitions in IRC +10216's envelope as observed 
with the IRAM 30 m telescope and Herschel/HIFI. The fine structure components are labelled with the respective /-transitions in 
the top of each panel, while the hyperfine components are indicated with vertical dotted lines. For the sake of clarity, we omitted the 
components that are too weak to be detected. See Sect. [4] for details on the spectroscopic structure of C2H. Additionally observed 
features related to other molecules are also identified. 



We refer to Table [2] for a summary of the presented data, and 
show an instructive comparison between the HIFI and SPIRE 
data of the C 2 H lines N = 6- 5uptojV = 9- 8in Fig. It 
is clear from these plots that the lower resolution of the SPIRE 
spectrum causes many line blends in the spectra of AGB stars. 



2.3. Herschel/PACS 

iDecin et alJ (l2010bh presented PACS data of IRC +10216, 
also obtained in the fram ework of the MESS programme 
dGroenewegen et aTl 1201 lb . The full data set consists of SED 
scans in the wavelength range 52 - 220 /an, obtained at differ- 
ent spatial pointings on November 12, 2009 (OD 182). We re- 
reduced this data set, taking into account not only the central 
spaxels of the detector, but all spaxels containing a contribution 
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Figure 2. SPIRE data in the range 475 - 795 GHz, containing the C2H-lines detected with HIFI. The shaded strips in the panels 
indicate the spectral ranges of the HIFI data shown in the panels of Fig. Q] The labels mark the principal contributors to the main 
spectral features. A comparison of emission in the range 609 - 6 14 GHz, as obtained with SPIRE and with HIFI, is shown in the 
inset panel. 
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Table 1. Summary of the IRAM 30 m and HIFI detections of C2H. Columns are, respectively, the instrument or band with which we 
detected the C2H emission, the main-beam efficiency 77MB, the half -power beam width HPBW, the integration time, the noise level 
in the spectra, the frequency resolution Av, the velocity resolution Av, the detected C2H rotational transitions, the frequency range 
in which we observed the lines, and the frequency-integrated intensity I* A = f T* k dv. 



Instrument 




HPBW 


Int.Time 


Noise 
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or Band 
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215 
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110.74 
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ncrsLtici/ nir 1 




















la 
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523941.9- 524000.4 
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lb 
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34.7 


13 
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0.25 
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611301.6-611364.1 


10.09 
9.44 


2a m 


0.75 
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698619.0 - 698738.0 


4.21 
4.47 


2b 
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15 


10 


6.0 


2.2 


N = 9- 8 
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785750.0 - 785875.0 <:!:) 
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Notes. ( , ) Only the horizontal polarisation was obtained. (:,:) Due to the low signal-to-noise ratio of this data set, these numbers are to be interpreted 
with caution. 



to the flux, i.e. 20 out of 25 spaxels in total. The here presented 
data set, therefore, reflects the total flux emitted by the observed 
regions, within the approximation that there is no loss between 
the spaxels. The estimated uncertainty on the line fluxes is of 
the order of 30%. Due to instrumental effects, only the range 
56 - 190//m of the PACS spectrum is usable. This range holds 
12 CO transitions J - 14 - 13 up to 7 = 42 - 41, covering energy 
levels from ~350cm _I up to ~3450crrr'. As is the case for the 
SPIRE data, line blends are present in the PACS spectrum due to 
the low spectral resolution of 0.08 - 0.7 GHz. 

2.4. IRAM 30 m: line surveys 

We combine the Herschel data with data obtained with the 
IRAM 30 m telescope a t Pico Veleta. The C 2 H N = 1 - tran- 
sition was observed by Kahane et alJ dl988l) : N — 2 - 1 was 
observed bv lCernicharo et al.l i20Mh N = 3- 2and/V = 4- 3 
were observed between January a nd April 2010, using the EMIR 
receivers as described in detail by dKahane et alll201ll in prep.). 
The data-reduction process of these different data sets is de- 
scribed in detail in the listed papers. The observed C2H tran- 



sitions are shown in Fig. [TJ and details on the observations are 
listed in TableQ] 

Table 2. Summarised information on the SPIRE spectrum con- 
taining the N = 7 - 6 transition of C2H. Av is the frequency reso- 
lution, Av is the velocity resolution, and cr is the rms noise. F v a* 
is the frequency-integrated flux in the given frequency range. 



Frequency range 


446 - 1575 GHz 


Integration time 


2664 s 


Av 


2.1GHz 


Av 


201 km s" 1 


F V ,A' 


16915 Jy MHz 


cr 


850 mJy 



3. The envelope model: dust & CO 

We assumed a distance d — 150pc to IRC +10216, in good cor- 
respondence with literature values (iGroenewegen et all 119981: 
iMen'shchikov et all 1200 it ISchoier et all 120071 and references 
therein). The second assumption in our models is that the effec- 
tive temperature T e $ = 2330 K, following the model of a large 
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set of m id-IR lines of C2H2 and HCN, presented bv lFonfrfa et alJ 
(120081) . We determined the luminosity L*, the dust-mass-loss 
rate Mdust, and the dust composition from a fit to the spec- 
tral energy distribution (SED; Sect. 13.11 ). The kinetic tempera- 
ture profile and the gas-mass-loss rate M were determined from 
a CO-line emission model (Sect. 13.2b . The obtained envelope 
model will serve as the basis for the CaH-modelling presented 
in Sect. [4] We point out that all modelling is performed in the 
radial dimension only, i.e. in ID, assuming spherical symmetry 
throughout the CSE. 

3.1. Dust 

The d ust modelling was performed using MCMax dMin et al.L 
l2009h . a Monte Carlo dust radiative transfer code. The best SED 
fit to the ISO SWS and LWS data, shown in Fig. [3] is based on a 
stellar luminosity L+ = \ 1300L(Ty This i s in goo d correspondence 
with the results o 

f|M 

en'shchikov et al. considering the 

difference in adopted distance. The combination of L* and 
gives a stellar radius of 20 milli-arcsec onds, in very good 
agreem ent with th e values reported by e.g. iRidewav & Keadvl 
(I19881 19 mas) and lMonnieret all (l2000t 22 mas). 

IRC +10216 is a Mira-t y pe pu lsator, with a pe- 
riod __of_6^9daxs_jLeBertti Il992h . Following Eq. 1 
of iMen'shchikov et all we find that L+ varies between 
i*,^=o~158OOL0 at maximum light (at phase ip-Q), and 
^*,^=o.5~6250Lo at minimum light (ip=0.5). Fig. [3] shows the 
SED-variability corresponding to the L* -variability. From the 
figure, it is visible that the spread on the photometric points can 
be accounted for by the models covering the full -range. 

The adopted dust composition is given in Table [3] The main 
constituents are amorphous carbon (aC), silicon carbide (SiC), 
and magnesium sulfide (MgS), with mass fractions of 53%, 25%, 
and 22%, respectively. The aC grai ns are assumed t o follow a 
distribution of hollow spheres (DHS; iMin et al 1 120031 with size 

0. 01 pm and a filling factor of 0.8. The population of the SiC and 
MgS dust grain s is represented by a continu ous distribution of el- 
lipsoids (CDE; iBohren & Huffman! 1 19831) . where the ellipsoids 
all have the volume of a sphere with radius «d=0. 1 pm. CDE and 
DHS are believed to give a more realistic approximation of the 
characteristics of circumstellar dust grains than a population of 
spherical grains (Mie-particles). Assuming DHS for the domi- 
nant aC grains was found to provide the best general shape of 
the SED. All dust species are assumed to be in thermal contact. 

The absorption and emission between 7 pm and 10pm 
and around 14pm that is not fitted in our SED model 
can b e explained by molecular bands of e.g. HCN and 
C2H2 dGonzalez- Alfonso & Cernicharol[T999l:ICernicharo et ail 
Il999h . 

Based on an average of the specific densities of the dust 
components p s = 2.41 gem -3 , we find a dust mass-loss rate 
of 4.0 x 10~ 8 Mq yr -1 , with the dust density pdust(/), dust tem- 
perature Td ust (r) and Qext/a^, with Q ext the total extinction effi- 
ciency, shown in Fig. |4] The inner radius of the dusty envelope, 

1. e. the dust condensation radius of the first dust species to be 
formed, is determined at Ri nna =2.7 R+ by taking into account 
pressure-dep endent condensatio n temperatures for the different 
dust species (iKama et all 120091) . These dust quantities are used 
as input for the gas radiative transfer, modelled in Sect. 13.21 



4 The extensive modelling of Men'shchi kov et al] was based on a 
light-curve analysis and SED modelling, and holds a quoted uncertainty 
on the luminosity of 20%. 



Table 3. Dust composition in the CSE of IRC +10216, used 
to produce the fit to the SED in Fig. [3] The columns list the 
dust species, the assumed shapes, the respective mass fractions, 
and the references to the optical constants of the different dust 
species. 



Dust species Shape Mass fraction References 

(%) 

Amorphous carbon DHS 53 j^reibischetaL Q993) 

Silicon carbide CDE 25 Pitman et al. (2008) 

Magnesium sulfide CDE 22 Beaemann et al. ( 1994) 




X (|im) 

Figure 3. SED models for IRC +10216, using M dust = 4.0 x 
10~ 8 MQyr~', and the dust composition and properties listed 
in Table [3] and spatial distribution and properties shown in 
Fig. g] The SED model is constrained by the ISO SWS and 
LWS data (ful l black); photometric poin ts (black diamonds ) 
are taken from iRamstedt et al.1 d2008l) and lLadial et al.l d2010l) . 
Dash-dotted grey: model at maximum light, with 15800 Lq, 
dashed grey: best fit to the ISO data, at phase ^=0.24, us- 
ing 1 1300L©, dash-triple-dotted grey: model at minimum light, 
with 6250 L©. The inset is identical, but plotted in log -log 
scale. 



3.2. CO 

To constrain the gas kinetic temperature T^ n (r) and the gas 
density p gas (r) throughout the envelope, we modelled the emis- 
sion of 20 rotational transitions of 12 CO (Sect. 13.21 ). measured 
with ground-based telescopes and Herschel/HIFI. The gas ra- 
diative transfer is treated with t he non-loca l thermal equ ilib- 
rium (NLTE) code GASTR0NO0M dDecin et all I2006L 1201 Pel) . To 
ensure consistency between the gas and dust radiative transfer 
models, we combine MCMax and GASTR0NO0M, by passing on the 
dust properties (e.g. density and opacities) from the model pre- 
sented in Sect. 13. II and Fig.|4]to the gas modelli ng. The general 
metho d behind this will be described in detail bv lLombaert et al.l 
(1201 2l in prep.). 

The large data set of high-spectral-resolution rotational tran- 
sitions of CO consists of 10 lines observed from the ground 
and 10 lines observed with HIFI, which are listed in Table [4] 
Since the calibration of ground-based data is at times uncertain 
dSkinner et all Il999h . large data sets of lines that are observed 
simultaneously and with the same telescope and/or instrument 
are of great value. Also, the observational uncertainties of the 
HIFI data are significantly lower than those of the data obtained 
with ground-based telescopes (20 - 40%). Since the HIFI data 
of CO make up half of the available high-resolution lines in our 
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Figure 4. Spatial distribution and properties of the dust for our 
best SED fit to the ISO data at <p=0.24 in Fig. [3j (a) dust den- 
sity PdustO), (b) dust temperature T dust (r) , and (c) Q ext /a d - The 
dashed lines indicate the inner and outer radii of the dusty enve- 
lope. 




10000 



Figure 5. Overview of the gas kinetic temperature and expan- 
sion velocity used in the radiative transfer model calculated with 
GASTRoNOoM. In the top panel, we indicated the exponents a 
from the T^ n (r) oc r"-power laws used to describe the kinetic 
temperature, and the radial ranges they apply to. See Sect. I3.2I 



To reproduce the CO lines within the observational un- 
certainties, we used L+ = 113OOL0, and a temperature profil^l 
that is a combination of three power laws: 7ki n (r) oc r -°- 5& 
for r < 9/?*, 7/^0) oc r 0A for 10/?* < r < 65 R+, and 
7\jnW oc r~ 12 at larger radii, ba sed on the work bv lFonfria et al.l 



(120081) and iDecinetalJ d2010al) . The minimum temperature in 
the envelope is set to 5 K. The radial profiles of T^ n (r) and u(r) 
are shown in Fig. [5] Using a fractional abundance CO/H2 of 
6 x 10~ 4 in the inner wind, we find that a gas-mass-loss rate 
M of 1.5 x lO~ 5 M0yr~' reproduces the CO lines very well. 
Combining the results from Sect. 13.11 with those from the CO 
model, we find a gas-to-dust-mass ratio of 375 , in the range of 
typica l values for AGB stars (10 2 - 10 3 ; e.g. iRamstedt et all 
12008b . 

The predicted CO line profiles are shown and compared 
to the observations in Fig. [6] All lines are reproduced very 
well in terms of integrated intensity, considering the re- 
spective observational uncertainties. The fraction of the pre- 
dicted and observed velocity-integrated main-beam intensities 
4iB,modei/4iB,data varies in the range 74 - 145% for the ground- 
based data set, and in the range 89 - 109% for the HIFI data 
set. The shapes of all observed lines are also well reproduced, 
with the exception of the IRAM lines (possibly due to excitation 
and/or resolution effects). 

Fig.Qshows the comparison between the observed CO tran- 
sitions / = 14 - 13 up to J - 42 - 41 in the PACS spectrum 
and the modelled line profiles, convolved to the PACS resolu- 
tion. Considering the uncertainty on the absolute data calibra- 
tion and the presence of line blends with e.g. HCN , C2H2, SiO, 
SiS, or H2O in several of the shown spectral ranges dDecin et all 
l20T0bh . the overall fit of the CO lines is very good, with excep- 
tion of the lines / = 14-13, 15- 14, and 16- 15. This significant 
difference between the model predictions and the data could, to 
date, not be explained. Both the HIFI data and the higher-/ PACS 
lin es are reproduced ver y well by the model. 

ISkinner et al.l JT999) noted the possible presence of a gradi- 
ent in the turbulent velocity u tul -b, with decreasing values for in- 
creasing radial distance from the star. We find no clear evidence 
in our data set that this decrease in i> tul -b is present. The model 
shown in Fig. [6] uses u tu ,b = 1.5 kms -1 , and reproduces the line 
shapes and intensities to a high degree. In Fig.|8]we compare the 
predictions of a set of CO-lines using different values for u tur b m 
order to assess the influence of this parameter. 
Assuming a lower value of 0.5kms _I or l.Okms -1 yields 
an incomplete reproduction of the emission profile between 
-44kms _1 and -41kms _I for all lines observed with HIFI. 
Assuming a higher value of 2kms~' leads to the overpredic- 
tion of the emission in this u-range for all observed lines, with 
the strongest effect for the lowest-/ lines. The influence of the 
different values of the turbulent velocity on the emission in the 
orange of -41 to -12kms _I is negligible. 



sample, these will serve as the starting point for the gas mod- 
elling. 

The adopted CO laboratory data are based on the work o f 
iGoorvitch & Chackerianl (| 1994) and IWinnewisser et al ] (ll997h . 
and are summarised in the Cologn e Datab ase for Molecular 
Spectroscopy (CDMS; iMuller et all 120051) . Rotational levels 
/ = up to / = 60 are taken into account for both the ground 
vibrational state and the first excit ed vibrational state. CO-H2 
collisional rates were adopted from lLarsson et al ] (l2002l) . 



4. C 2 H 

4. 1 . Radical spectroscopy 

C2H (• C = C — H) is a linear molecule with an open shell con- 
figuration, i.e. it is a radical, with a n electronic 2 £ + ground-state 
configuration (Miiller et al., 2000). Spin-orbit coupling causes 
fine structure (FS), while electron-nucleus interaction results 
in hyperfine structure (HFS). Therefore, to fully describe C2H 
in its vibrational ground-state, we need the following coupling 

5 The central star is assumed to be a black body at a temperature T s g . 
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Figure 6. Comparison of the observed 12 CO emission lines (black histogram), and the line profiles predicted by the GASTRoNOoM- 
model (red dashed line) with parameters as listed in Table [5] The line transitions J - (J - 1) and the telescopes with which they 
were observed are indicated in the upper left corner of every panel. The factor /MB.modei/^MB.data is given in the upper right corner of 
every panel. 



scheme: 

J = N + S (1) 
F - J + I, (2) 

where N is the rotational angular momentum, not including the 
electron or quadrupolar angular momentum (S and /, respec- 
tively), J is the total rotational angular momentum, and F is the 
nuclear spin angular momentum. 

The strong AJ = AN FS components have been detected 
up to N — 9 - 8. Relative to the strong components, the weaker 
AJ = components decrease rapidly in intensity with increasing 
N, and they have been detected up to N — 4 - 3. The HFS has 
also been (partially) resolved for transitions up to N = 4 - 3. 

Spectroscopic properties of the ground vibr ational state of 
C2H h ave most recently been determined by iPadovani et al.l 
(120091) and are the basis for the entry in the CDMS. 

Our treatment of the radiative transfer (Sect. 14.3b does not 
deal with the FS and HFS of C2H, and is limited to the predic- 
tion of rotational lines N — » N' with AN = 1, according to a S 
approximation. Therefore, all levels treated are described with 
only one quantum number N. The final line profiles are calcu- 
lated by splitting the total predicted intensity of the line over 
the different (hyper)fine components, depending on the relative 
strength of these components (CDMS), and preserving the total 
intensity. Note that, strictly speaking, this splitting is only valid 
under LTE conditions, but that, due to our approach of C2H as a 
'Z-molecule, we will apply this scheme throughout this paper. 



C2H has one bending mode vo, and two stretching modes, 
v\ and V3. The rotational ground-state of the bending mode (V2) 
is situated ~530K above the ground state. The C — C stretch 
(v 3 ) at -2650 K is strong. The C - H stretch (v,) at -4700 K is 
weak, and is resonant with the first excited 2 II electronic state at 
~5750 K of C2H. For each of these vibrational states — ground- 
state, V2, V3, and v\ — we consider 20 rotational levels, i.e. rang- 
ing from N = up to N = 19. 

The equilibrium dipole moments of the ground and 
first excited electronic states have been calcu la ted as /u = 
0.769 and 3.004 Debye, respectively (IWoonL Il995h . The 
dipole moments of the fundamental vibrational states in the 
ground electronic state are also assumed to be 0.769 Debye. 

This is usually a good assumption as shown in the 

case of HCN by iDeleon & Muentej dl984h . The band 
dipole moments for rovibrational transiti ons were calcu- 
lated from infrared intensities published by iTarroni & Carted 
(12004 : //(v 2 = l->0) = 0.1 10 Debye, ju(v 3 = l->0) = 0.1 78 Debye, 
and /4vi = l— >0) = 0.050 Debye. The influence of the vibra- 
tionally excited states on transitions in the vibrational ground 
state will be discussed in Sect. 14.31 

A more complete treatment of C2H will likely have to take 
into account, e.g., overtones of the v?-state. These h ave non- 
negligible intrinsic strengths ( ITarroni & CarteH 120041) because 
of anharmonicity and vibronic coupling with the first excited 
electronic state. Spectroscopic data for several of these are al- 
ready available in the CDMS. A multitude of high-lying states 
may also have to be considered as they have rather large intrin- 
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Figure 7. Comparison of CO-emission lines measured with PACS (black histogram) and the predicted line profiles (red dotted line). 
The transitions / - (/ - 1) are labelled in the top left corner of each subpanel. 



sic intensities because of the vibronic interaction between the 4.2. Observational diagnostics 
ground and the first excited electronic states. 

The C2H-emission doublets due to the molecule's fine struc- 
ture (Sect. 14. U are clearly present in all high-resolution spectra 
shown in FigQ] For = 1 -0, 2—1, and 3-2, the hyperfine struc- 
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Figure 8. Zoom on the blue wing of a selection of the CO line profiles, with the transitions J -(J - 1) labelled in the top left corner 
of each panel. We show the comparison between the observed lines and predicted lines with assumed values of u tu ,b =0.5kms _1 
(full grey), 1.0 km s -1 (dotted grey), the adopted value of 1.5 km s (dashed grey), and2.0kms~' (dash-dotted grey). 



ture is clearly detected in our observations. The double-peaked 
line profiles are ty pical for spatially resolved o ptically thin emis- 
sion (Olofsson, in lHabing & Olofssonl 12001 . and indicate that 
the emitting material has reached the full expansion velocitjfl 
i.e. Vco = 14.5 km s~'. This is in accor dance with the obse rved 
position of the emitting shells (~16"; iGuelin et al.L 1 1999b and 
the velocity profile derived for this envelope (Sect. l3~2l i. 

To derive information on the regime in which the lines 
are ex cited, we construct a rotational diagram dSchloerb et al.L 
119831) . using 



N 



Z(T rot ) 



x exp 



3k x 10 36 7„, colT 
8tt 3 vSfi 2 



(3) 



where TV is the column density, Z the temperature dependent par- 
tition function, E u the energy of the upper level of the transition 
in cirT 1 , k the Boltzmann constant in units of cirr 1 Kr 1 , T rot the 
rotational temperature in K, /i the dipole moment in Debye, and 



v the frequency of the transitions in Hz. 
integrated antenna temperature, 



/bff^MB 



dv 



is the velocity- 



(4) 



corrected for the beam filling factor (lKrameit[T997l) 

_ / 1 — exp |— |^ Vln 2 X $source/$beam]"j » ^source ^ ^beam 
\ 1 > ^source — ^beam 

and for the beam efficiency ^mb (see TableQ}. We have assumed 
a uniform emission source with a size of 32" in diameter for all 
C2H transitions, based on the interferometr ic observations of the 
N = 1 — emission bv lGuelin et all (fl999l) . 

From Fig. [9] we see that the C2H emission is characterised 
by a unique rotational temperature 7^ = 23.3(+1.3)K, and a 
source-averaged column density TV = 3.84(±0.09) x 10 15 cm -2 . 



6 We assume that the half width at zero level of the CO emission 
lines, i.e. 14.5 km s~' , is indicative for the highest velocities reached by 
the gas particles in the CSE. 



N = 3.84 (0.09) x 10*' 5 cnr* 
T rot = 23.4(1 .3) K 




50 100 

EJk (K) 



150 



Figure 9. Rotational diagram of the measured C2H transitions, 
based on Eq. [3] Transitions — in the 'Z approximation — are 
labelled in grey. The numbers given in parentheses are the un- 
certainties on the derived column density and rotational temper- 
ature. 



The unique rotational temperature points to excitation of the 
lines in one single regime, in accordance to the abundance peak 
in a c onfined area that can be characterised with one gas temper- 
ature dGuelinet aUll999l) . 



4.3. Radiative transfer modelling 
4.3.1. Collisional rates 

For rotational transitions within the vibrational ground state and 
within the — 1 bending state and the v\ — 1 and V3 = 1 
stretching states we adopted the recently pu blished collisional 
rates of HCN-H 2 bv lDumouchel et all (l2010h . The lack of colli- 
sional rates for C2H, and the similar molecular mass and size of 
HCN and C2H inspire this assumption. The collisional rates are 
given for 25 different temperatures between 5 and 500 K, hence 
covering the range of temperatures relevant for the C2H around 
IRC +10216. Since v 3 and v { are at -2650K and -4700 K, re- 
spectively, collisional pumping to these vibrationally excited 
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Table 4. Overview of the transitions of CO and C2H 
shown throughout this paper. We list the transitions, the tele- 
scope/instrument, the dates of observation, and the litera- 
ture references in case of previo u sly pu blish ed data. These 
referen ces are: (1) Huggins et alj ( 1988b. (2) Olpfsson et alj 
d!993l), (3) ICernicharo J. and Kahane et all j201lL in p repX 



(41 iGroenewegen et all (Il996l) . (5 ) IWang et all Jl994l) . (6) 
iKahaneetalJ d20lU in prep.), (7 ) ITevssier et al] d2006l). (8 ) 
this paper. (9) iDecin et alj d20 10b h (10") IKahane et alJ d!98 



(1 1) ICernicharo et all d200oh . (121 lAgundez et al.1 d201( 
ICernicharo et al. ( 201 lal) . 



Transitions 
,2 CO(J -> J- 1) 

1-0 
1-0 
1-0 

1- 

2- 1 

3- 2 
3-2 

3- 2 

4- 3 
6-5 



Telescope 



Date 



5-4 ... 


11 - 


10 


14-13 ... 


16- 


15 


14-13 ... 


25- 


24 


27-26 ... 


42- 


41 


C 2 H (N —> N 


- 1) 





1- 

2- 1 

3- 2 

4- 3 

6- 5 

7- 6 

8- 7 



Ref. 



NRAO 


Jun, 1986 


1 


SEST 


Oct 13, 1987 


2 


IRAM 


Sept 3, 2004 


3 


IRAM 


Oct 8, 1991 


4 


IRAM 


Aug 1, 2003 


3 


CSO 


Jun, 1993 


5 


JCMT 


Jul 13, 1992 


4 


IRAM 


Feb 3, 2010 


6 


CSO 


Feb 2002 (a) 


7 


CSO 


Feb 2002 (fl) 


7 


HIFI 


May 11-13, 2010 


8 


PACS 


Nov 12, 2009 


8,9 


IRAM 


May 27, 1985 


10 


IRAM 


Jan 7, 2002 


11 


IRAM 


Jan 29, 2010 


12 


IRAM 


Mar 28, 2010 


6, 13 


HIFI 


May 11-13, 2010 


8 



Notes. <n) ITevssier et al.l d2006t) state that observations were carried out 
between September 2001 and February 2002. A search of the CSO 
archive database resulted in observations in February 2002 that could 
be linked to their paper. 



Table 6. Properties of the transitions N — > N' used to model 
the C2H lines presented in Fig. Q] v is the rest frequency of the 
transition, E^jk is the energy level of the upper level N of the 
transition divided by the Boltzmann constant k, A^^n' is the 
Einstein-A coefficient of the transition, S n^n' is the theoretical 
line strength of the line, is the statistical weight of the upper 
level N. 



N 


N' 


V 

(MHz) 


E N /k 
(K) 


Aat-»w 
(s" 1 ) 




5 /v_,/v< 


gN 


1 





87348.635 


4.2 


1.53x10- 


-6 


1 


3 


2 


1 


174694.734 


12.6 


1.47x10- 


-5 


2 


5 


3 


2 


262035.760 


25.2 


5.31x10- 


-5 


3 


7 


4 


3 


349369.178 


41.9 


1.30x10 


-4 


4 


9 
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5 


524003.049 


88.0 


4.57x10- 


-4 
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13 
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611298.435 


117.4 
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Figure 10. Fractional abundance of 12 CO (black) and C2H (red) 
relative to H2: model including the density enhancements (full 
lines), and models without den sity enhancem ents, assuming the 
average interstellar UV field of lDraind (1 19781) scaled by a factor 
1 (dotted lines), 0.5 (dashed lines), and 0.25 (dash-dotted lines). 
See Sect. I4.3.2l for a discussion of this plot. 



Table 5. Parameter set used to model the CO lines, as discussed 
in Sect. [12] 



d 


150pc 


M dusl (n) 


4.0 x 10~ 8 M Q yr" 1 




113OOL 


M gas 


1.5 x lO^Moyr 1 


R+ 


4.55 x 10 13 cm 


c6/h 2 


6x 10- 4 




* 655 R Q 




14.5 km s- 1 


R <«) 

Jv inner 


2.7 R, 


»turb 


1.5 kms 1 


^onter 


25000 tf* 


7kin(r) 


oc r-° 58 (ltf* <r< 9R t ) 




2330 K 




oc r-° 40 (9tf, < r< 65 R+) 








oc r - 120 (65^ < r) 



Notes. {fl) From SED modelling. 



levels is unlikely for the low temperatures prevalent in the ex- 
citation region of C2H, i.e. the radical shell at ~16". Hence, 
for rovibrational transitions between the vibrationally excitated 
states (V2 = 1, V3 = 1, v\ — 1) and the vibrational ground state 
we adopted the same rates as for collisionally excited transi- 



tions within the vibrational ground state, but scaled down by a 
factor 10 4 , comparable to what ha s been done e.g. f or H2 O by 
Deguchi & Nguven-O-Rieul (Il990h : iMaercker efall d2008l) and 
Decin et al.l d2010a) . We did not consider collisionally excited 
transitions between the excited vibrational states. 

We note that recent calculations bv lDumouchel et al] (1201 Ol) 
have shown that the HNC collision rates at low temperatures dif- 
fer by factors of a few from those of HCN. Since we may ex- 
pect similar errors for C2H, there is the need for accurate C2H 
— H2 collisional rates as well. The gas density, therefore, is not 
constrained better than this factor in the models presented in 
Sects. Ll2land l4T2l 



4.3.2. Constraining the C2H abundance profile 

The calculation of the C2H/H2 fr actional abundance is based on 
the chemical model discussed bv lAgundez et al. Iipoioh assum- 



ing the average interstellar UV field from lDraine dl978h . and a 
smooth envelope structure. The CO and C2H abundances corre- 
sponding to the envelope model presented in Sect. l3.2l are shown 
in Fig. [TO] 
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The Pla t eau d e Bure Interferometer (PdBI) maps of 
IGuelin et al.l (1 19991) show that the 3mm-lines of the three rad- 
icals C2H, C4H, and CeH have their brightness peaks at a ra- 
dial angular distance of ~16" fro m the central s tar. In Fig. 1 1 l(a)| 
the contours of the C2H map by IGuelin et al.l d 19991) are over- 
layed on the normalised brightness distribution of the N = 1 — 
transition predicted by a (one-dimensional) LTE model based on 
the abundance profile mentioned above. The extracted PdBI con- 
tours correspond to the velocity channel at the systemic velocity 
of IRC +10216, representing the brightness distribution of the 
CiH-transition in the plane of the sky. From Fig. 1 1 l(a)| it is clear 
that this model leads to a brightness distribution that is concen- 
trated too far inwards to be in agreement with the shown PdBI 
contours. 

Several authors have established that the mass-loss process 
of IRC +1 0216 shows a complex time -dep endent beha v iour. T he 
images of iMauron & Hugginsl d 19991) and iLeao et all d2006l) in 
dust-scattered stellar and ambien t (optical) light, and the CO(J = 
1-0) maps of lFong et alJ (l2003h show enhancements in the dust 
and the gas densi ty, with quasi-periodic behaviour. Recently, 
iDecin et all d2011h showed the presence of enhancements out to 
~ 320" based on PACS-photo metry, consistent with the images 
of lMauron & Hugginsl dl999l) . 

iBrown & Millarl d2003l) presented a chemical envelope 
m odel incorporating the dus t density enhancements ob served 
bv IMauron & Hugginsl dl999t) . ICordiner & Millarl d2009l) added 
density enhancements in the gas at the positions of the dust den- 
sity enhancements. Th ey assumed complete dust-ga s coupling, 
based on the work by iDinh-V-Trung & Liml d2008l) who com- 
pared map s of the molecular sh e lls of HC3N and HC5N with the 
images bv IMauron & Hugginsl dl999l) . To mimic these density 
enhancements in our model and evaluate the impact on the exci- 
tation and emission distribution of C2H, we use the approxima- 
tion of ICordiner & Miliar! d2009l) . adding (to our "basic" model, 
presented in Sect. 13. 2t ten shells of 2" width, with an intershell 
spacing of 12", starting at 14" from the central star. These shells 
— located at 14"— 16", 28"-30", etc. — are assumed to have 
been formed at gas-mass-loss rates of a factor six times the rate 
in the regions of "normal density", the intershell regions. We 
note that similar episodic mass loss was not added to the "ba- 
sic" dust model presented in Sect 13.11 The uncertainties on the 
data presented in Sect. 13.11 and the photometric variability of 
IRC +10216 are too large to constrain the expected small effect 
of the inclusion of these enhancements. 

These density enhanc e ments are included in the chemical 
model of lAgundez et alJ d2010h by combining two regimes. 
Firstly, a model is run in which the chemical composition of 
a parcel of gas is followed as it expands in the envelope, con- 
sidering an augmented extinction (Ay) contribution from the 
density-enhanced shells located in the outer CSE. The density 
of the parcel is determined by the non-enhanced mass-loss rate 
(M). A second model is run in which the composition of a 
density-enhanced shell is followed as it expands. The final abun- 
dance profile follows the abundance from the first model for the 
non-enhanced regions and follows the abundance of the second 
model for the density-enhanced shells. The resulting abundance 
profiles for 12 CO and C2H are compared to those correspond- 
ing to a smooth model without enhancements in Fig. [TO] The 
12 CO-abundance is affected by the inclusion of density enhance- 
ments only in the outer regions of the CSE. The impact on the 
radiative-transfer results for 12 CO is negligible, with differences 
in modelled integrated intensities of at most 3% compared to the 
model from Sect. 13.21 



We find that the presence of density enhancements signifi- 
ca ntly alters the abundance p rofile of C2H, as was also discussed 
bv lCordiner & Millarl d2009l) . The abundance peak at -500 in 
the smooth model, is shifted to ~800 in the enhanced model, 
as shown in Fig. [10] This shift indicates that the photochemistry, 
in this particular case the photodissociation of C2H2, is initiated 
at larger radii than in the smooth model. This is caused by the 
stronger extinction due to the density enhancements in the outer 
envelope. At the position of the two innermost density enhance- 
ments the fractional C2H abundance is lower than in the inter- 
shell regions. The absolute abundance of C2H in these regions, 
however, is higher by a factor 2-3 than in the neighbouring in- 
tershell regions. This is a combined effect of (1) the augmented 
shielding of C2H2 from incident interstellar UV radiation (which 
is no longer true for the outer shells), and (2) chemical effects 
such as faster reactions with e.g. C2H2 to form larger polyynes. 

We note that the correspondence between the interferomet- 
ric observations and the modelled number densities has signifi- 
cantly improved by including density enhancements (Fig. |l 1(b)) . 
However, the peak intensity is located at somewhat too small 
radii compared to the observations. The adopted enhancements 
in the envelope are linked to M-values of a factor six times as 
high as for the intershell r egions, in correspondence with the val- 
ues used in the model of ICordiner & Millarl d2009l) . Increasing 
this factor or the number of shells in our model would increase 
the total mass of material shielding the C2H2 molecules from 
photodissociation by photons from the ambient UV field, and 
would hence shift the peak intensity of the predicted emission 
outwards. Furthermore, the presence of nu merous dust arcs out 
to ~320" is reported bv lDecin et al. I d20ll and their effect will 
be included in future chemical models. 

The "basic" model from Sects. 13.21 and 14.31 and the above 
introduced model wit h the density en hancements assume the in- 
terstellar UV field of iDraind dl978l To assess the effect of a 
weaker UV field on the photodissociation of C2H2, and hence on 
the spatial extent of C2H, we tested models assuming interstellar 
UV fields weaker by a factor of 2 and 4. The corresponding C2H 
abundance profiles are shown in Fig. [TUJ and the the modelled 
N = 1 — brightness distribution for the latter case is shown 
in Fig. |1 1(c)] Although the correspondence in this figure is very 
good, the C4H and CeH abundances produced by this chemi- 
cal model cannot account for the c oinciding brightness peaks of 
C 2 H, C 4 H, and C 6 H as reported bv IGuelin etaiT d 1 9991) . In con- 
trast, the density-enhanced model does reproduce this cospatial 
effect. 

4.3.3. Excitation analysis 

Assuming the C2H abundance obtained from the envelope model 
with density enhancements, we tested the influence of includ- 
ing the vibrational modes by modelling four cases: (1 ) including 
only the ground-state (GS) level, (2) including GS and V2 levels, 
(3) including GS, V2, and V3 levels, and (4) including GS, V2, V3, 
and vi levels. An overview of the ratio of the predicted and the 
observed integrated intensities, /Modei/^Data, for these four cases 
is shown in Fig. [T2] The comparison of the line predictions for 
case (1 ) and case (4) under NLTE conditions is shown in Fig.fOI 
Under LTE conditions, /Model /feata for a given transition is 
the same for all four cases, indicating that the vibrationally ex- 
cited states are not populated under the prevailing gas kinetic 
temperatures of ~20 K in the radical shell. Under NLTE condi- 
tions, however, the predicted line intensities are very sensitive 
to the inclusion of the vibrationally excited states. In the two 
cases where the V3-state is included, the predicted intensity of 
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interstellar UV field a quarter as strong as in the standard model 
presented in panel (a) 



Figure 11. Contours (in white) of the brightness distribution of the C2H N = 1 — transition measured by iGuelin et al.l (1 19991) 
with PdBI, overlayed on the normalised predicted brightness distribution, assuming LTE conditions. The three panels represent the 
results obta ined with three different chemical models. For the map in panel (a), we assumed a constant mass loss, and the interstellar 
UV field of lDraind(fl978l).The ma p in pa nel (b) shows the brightness distribution when assuming density enhancements in the CSE 
as described bv lCordiner & Millarf (120091) and as discussed in Sect. 14.3.2] For panel ( c), we assumed again a constant mass loss, but 
an interstellar UV field that is only 25% as strong as that presented by Draind(ll978l) . 



the N = 1 — transition is only 9% of the observed value, while 
transitions with N up > 4 are more easily excited. This is linked to 
the involved Einstein-A coefficients. For example, Table|7]gives 
the Einstein-A coefficients of the transitions involving the N - 2 
level of the vibrationally excited states and the N - 3 and N = 1 
levels of the ground state. The higher values of the Einstein-A 
coefficients for transitions to the N - 3 level lead to a more 
effective population of this level than of the N - I level. In par- 
ticular, when including V3 and v\ , this causes an underprediction 
of the intensity of the N - \-Q ground state transition, while this 



effect is insignificant when only V2 is included. A visualisation 
of this scheme is shown in Fig. [14] This pumping mechanism is 
not limited to AT = 3, but also affects higher levels, explaining 
the easier excitation of the transitions with N up > 4 mentioned 
before. 



4.3.4. Turbulent velocity 

As was done for CO, we tested the influence of y t mb on the 
predicted C2H-emission, considering values u tur b=0.5, 1.0, 1.5, 
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Figure 13. Comparison of the measured C2H spectra (black histograms; see Sect, \2\and Fig. \I\for a description and for the identifi- 
cation of additional spectral features) and GASTR0NO0M model predictions under NLTE conditions for the case where (red) only 
the ground state is included, and (blue) the ground state, and the three vibrational modes are taken into account. These predictions 
are based on the "enhanced" abundance profile shown in Fig.[10]and the 'E-approximation of C2H. The transition N - (N - 1) and 
the ratio 7 mo dei/^data are stated in the upper right corner of each panel, following to the colour code of the plots. 



and 2.0kms~'. We find that the predicted line intensity in- 
crease^] with 25-55% when v ta ± is increased from 0.5 km s _I to 
2.0kms~ 1 . From the comparison of the predicted and observed 
shapes of the emission lines, we conclude that values of u tur b in 
the range 0.5 - 1.5kms _1 give the best results. This is consistent 

7 This is valid for all lines, except for the N = 1 - transition, where 
the predicted intensity decreases by 40%. 



with our discussion in Sect. l3.2l for t he CO-lines, and wit h values 
generally used for AGB envelopes dSkinner et alj. 1 1999b . 

4.4. Vibrationally excited states 

The removal of a hydrogen atom from C2H2 through photodis- 
sociation causes the produced C2H molecule to be bent, rather 
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Figure 12. Overview of the ratio /Model /^Data of integrated in- 
tensities for the C2H lines in Fig. Q] The x-axis is labelled ac- 
cording to the upper AMevel of the transition, N up . Different 
approximations of the C2H molecule are plotted with different 
symbols: crosses (+) for the model including only the ground 
state (GS), asterisks (*) for the model including GS and the V2 
state, diamonds (o) for the model including GS, V2, and V3, and 
triangles (A) for the model including the GS and all three vibra- 
tionally excited states vn, V3, and v\. Black symbols connected 
with full lines, and red symbols connected with dotted lines rep- 
resent NLTE and LTE models, respectively. 

Table 7. Einstein-A coefficients (in s~') for the transitions in- 
volving the N - \ and N - 3 levels of the ground state (GS). 









lower level 








GS, N = 3 


GS, N =1 
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8.16 x 10~ 2 


upper level V3, 
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= 2 


3.67 x 10 1 


2.50 x 10 1 


Vl, 
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= 2 


1.68 x 10 1 


1.13 X 10 1 



than linear dMordaunt et al.L [l998h . This means that this forma- 
tion route of C2H favours the population of the V2 = 1 state over 
the population of the ground state. However, the Einstein-A co- 
efficients for rovibrational transitions in the band V2 = 1 — > 
are in the range 0.01 - 1 s _1 , and by far exceed the photodis- 
sociation rates of C7H7, which are of the order of 10~ 9 s _1 
(I van Dishoeck et aTll2006h . Hence, we did not take this V2-state 
population effect into account. 

Recently, iTenenbaum et al.l d2010b reported on the detection 
of the Coli N = 3-2 transition in the vn — 1 state, with 
the Arizona Radio Observatory Submillimeter Telescope (ARO- 
SMT). Observed peak intensities are of the order of ~5 - 10 mK 
in antenna temperature, with reported noise levels around 3 mK, 
and a total velocity-integrated intensity of ~0.2Kkms _1 . The 
V2 = 1, N = 3 level is located at ~560K and is not populated 
under LTE conditions, given the gas-kinetic temperature in the 
radical shell (~20K). Under NLTE conditions, however, the v%- 
levels are easily populated. The prediction of the V2 = 1, N = 
3-2 emission under the different 'E-approximations is shown 
in Fig. [15] Under the 1 ^-approximation, and assuming an HPBW 
of 29" at 261 GHz for ARO, we predict a U-shaped profile with 
peak antenna temperature ~3 mK. Taking into account the ob- 
served fine structure, we hence predict an integrated intensity 
which is a factor ~4 lower than the observed value. Considering 
the substantial uncertainties (the low signal-to-noise ratio of the 
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N=2 
N= 1 
N= 




GS 



Figure 14. Visualisation of the transitions responsible for the de- 
population of the N = 1 level in the vibrational ground state 
(GS), where vibrational states are indicated on the horizontal 
axis, and rotational levels on the vertical axis. The colour of 
the arrows representing the transitions is indicative of the mag- 
nitude of the Einstein-A coefficients of these transitions, as given 
by the colour bar. 



GS 
GS 



GS + v 2 + v 3 + v, 




Figure 15. The predicted line profile for the N - 3 - 2 transition 
in the v% = 1 vibrationally excited mode, under assumption of 
NLTE conditions, and different 'E-approximations, as it would 
be observed with the Arizona Radio Observatory Submillimeter 
Telescope. 



ARO observations, and the limited C2H approximation in our 
model), this agreement is satisfactory. 



5. Summary 

We presented new data of CO and C2H obtained with HIFI, 
PACS, SPIRE and the IRAM 30 m telescope. High-resolution 
spectra of CO transitions up to / = 16-15, and of C2H transi- 
tions up to = 9 - 8 were presented and modelled. The HIFI 
data of both CO and C2H are the first high-frequency -resolution 
detections of these lines. They were obtained in the framework 
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of a spectral survey of IRC +1 0216 over the complete frequenc y 
range of the HIFI instrument (ICernicharo et all 1201 Obi 1201 lbl) . 

From an SED fit to ISO data, PACS data, and a set of pho- 
tometric points, covering the wavelength range 0.1 - 1000 /im, 
we obtained a dust-mass-loss rate of 4.0 x 10~ 8 M0yr~' and a 
luminosity of 11300L© at a distance of 150pc. This luminos- 
ity value is valid at tp = 0.24, the phase at which the ISO data 
were obtained. The luminosity is then expected to vary between 
6250 Lq and 15800 L© throughout the star's pulsational cycle, 
which has a period ~649 days. 

In order to model IRC+10216's wind, we performed the 
radiative transfer of the dusty component of the CSE consis- 
tently with the gas-radiative transfer of CO. The set of 20 high- 
spectral-resolution CO lines was modelled in order to constrain 
the physical parameters of IRC +10216's CSE. The kinetic tem- 
perature through out the envelope was desc ribed by previous re - 
sults reported by iFonfria eTai] d2008l) and lDecin etal] d2010al) . 
and is now further constrained by the combination of all the 
high-resolution CO lines. The temperature profile is charac- 
terised by r kin (r) oc r-°- 58 for r < 97?*, T kin (r) oc r 0A for 
10 R+ < r < 65/?*, and T^ D (r) oc r~ 12 at larger radii, with an 
effective temperature 7^=2330 K. The derived mass-loss rate 
is 1. 5 x 10~ 5 Mr?) yr" 1 . This is consistent with earlier results 
(e.g. ICernicharo et all 120001) and gives a gas-to-dust-mass ra- 
tio of 375, in line with typical values stated for AGB stars (e.g. 
iRamstedt et all 120081) . Furthermore, we showed a very good 
agreement between the predictions for CO lines up to J = 42-41 
and the newly calibrated PACS spectrum of IRC + 10216. It is the 
first time that such a large coverage of rotational transitions of 
CO is modelled with this level of detail. 

We extended our envelope model wit h the inclusion of 
episod ic mass loss, based on the model of ICordiner & Miliar! 
(120091) . This assumption proved very useful in reconciling the 
modelled C2H emission with t he PdBI map of the iV = 1 - 
transition of C2H published by iGuelin et al.1 dl999l) . and in re- 
producing the observed line intensities. A decrease of a factor 
four in the strength of the interstellar UV field also leads to a 
satisfactory reproduction the PdBI map, but resulted in poorly 
modelled line intensities. The inclusion of density enhancements 
in IRC+10216's CSE is further supported by observational re- 
sults based on maps of dust -scattered light (|Mauron & HugginsL 
[1991 . molecular emisson dFong et ail, 12003b, and ph otometric 
maps recently obtained with PACS dDecin et allfeoi lb . 

The ground-based observations of C2H transitions involve 
rotational levels up to N = 4 with energies up to ~17 cm -1 , cor- 
responding to temperatures ~25 K. This temperature is close to 
the gas kinetic temperature at the position of the radical shell 
(~20K). The recent detection of strong C2H-emission involving 
levels with energies up to ~150K, however, calls for an efficient 
pumping mechanism to these higher levels. Due to the spectro- 
scopic complexity of C2H, with the presence of fine structure 
and hyperfine structure, we approximated the molecule as a 1 S- 
molecule, exhibiting pure rotational lines, without splitting. We 
illustrated that the inclusion of the bending and stretching modes 
of C2H is crucial in the model calculations, since high-energy 
levels are much more efficiently (radiatively) populated in this 
case. At this point, we have not yet included overtones of the 
vibrational states, nor did we treat the resonance between vi- 
brational levels in the electronic ground state and the first elec- 
tronically excited A 2 LI-state. Applying our simplified molecu- 
lar treatment of C2H, we can explain the strong intensities of 
the rotational lines in the vibrational ground state, except for the 
N = 1 — transition. We are also able to account for the excita- 



tion of the recently observed rotational transition in the — 1 
state, showing the strength of our approach. 
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